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Abstract 

The paper revises the explicit integration of the classical Steklov-Lyapunov systems via sep- 
aration of variables, which was first made by F. Kotter in 1900, but was not well understood 
until recently. We give a geometric interpretation of the separating variables and then, apply- 
ing the Weierstrass hyperelliptic root functions, obtain explicit theta-function solution to the 
problem. We also analyze the structure of its poles on the corresponding Abelian variety. This 
enables us to obtain a solution for an alternative set of phase variables of the systems that has 
a specific compact form. 



1 Introduction 

The motion of a rigid body in the ideal incompressible fluid is described by the classical Kirchhoff 
equations 

r > r ^ dH dH . dH 

K = K x Ti? +P X TT' P = P x u-> 
UK op op 

where K, p € R 3 are the vectors of the impulsive momentum and the impulsive force, and H = 
H(K,p) is the Hamiltonian, which is quadratic in K, p. Note that this system always possesses two 
trivial integrals (Casimir functions of the coalgebra e*(3)) (K,p), (p,p) and the Hamiltonian itself is 
also a first integral. 

Steklov [2U] noticed that the classical Kirchhoff equations are integrable under certain conditions 
i.e., when the Hamiltonian has the form 



1 3 

Hl = 2 E (^^ + 2 ^^7 A ^P« + " 26 «(^ - 6 7) 2 P«) > (a,/3, 7 ) = (1,2,3), (1) 



"AMS Subject Classification 37J35, 70E40, 70H06, 70G55, 14H40, 14H42 



1 



bi,b2,bs and v being arbitrary parameters. Under the Steklov condition, the equations possess 
fourth additional integral 



1 3 

4 = r^fc- 2vb a K aVa + v\bp - b^fpl) . (2) 



3 

2 

a = l 

Later Lyapunov [17] discovered an integrable case of the Kirchhoff equations whose Hamiltonian 
was a linear combination of the additional integral ((2]) and the two trivial integrals. Thus, the 
Steklov and Lyapunov integrable systems actually define different trajectories on the same invariant 
manifolds, two-dimensional tori. This fact was first noticed in [14) . 

In the sequel, without loss of generality, we assume v = 1 (this can always be made by an 
appropriate rescaling p — > p/v). 

The Kirchhoff equations with the Hamiltonians |T]) , ([2]) were first solved explicitly by Kotter [16] , 
who used the change of variables (K,p) — > (z,p): 

2z a = K a - (bp + 6 7 )p Q , a = 1,2,3, (a, /3, 7) = (1,2,3), (3) 

which transforms the Steklov-Lyapunov systems to the form 

z = z x Bz — Bp x Bz , p = pxBz, B = diag (61, 62, 63) (4) 

and, respectively, 

z = p x Bz , p — p x (z — Bp) . (5) 

Kotter implicitly showed that the above systems admit the following Lax representation with 
3x3 skew-symmetric matrices and a spectral parameter 

L(s) = [L(s),A(s)], L(s),A(s) € so(3), sGC, 

(6) 

where £ Q/ 3 7 is the Levi-Civita tensor. Equations Q and © are generated by the operators 

A{s) al 3 = '^ 1 \J (s - b a )(s - bp) bjZj , resp. A(s) a p = e a pjyj (s - b a )(s - bp)p 1 . (7) 

The radicals in ©-([7]) are single-valued functions on the elliptic curve £ , the 4-sheeted unramified 
covering of the plane curve £ — {w 2 = (s—bi)(s—b2)(s—bs)}. For this reason, the Lax representation 
has an elliptic spectral parameter. 

Writing out the characteristic equation for L(s), we arrive at the following family of quadratic 
integrals 

3 

F(s) = ^(s - & 7 )(z 7 + s Pl ) 2 = Jis 3 + J 2 s 2 + 2sH 2 - 2Hi , (8) 
7=1 

where 

H! = ±{z,Bz), H 2 = ±{z,z)-{Bz,p), J 2 =2{z,p)-{Bp,p), Ji = (p,p). (9) 

It is seen that under the Kotter substitution ([3]) the functions J\ , J2 transform into invariants 
of the coalgebra e*(3), whereas the integrals H\(z,p), H%(z,p) (up to a linear combination of the 
invariants) become the Hamiltonians ([T]),([2]). 

An analog of the Lax pair ([6]) was later rediscovered in and was used to obtain theta-function 
solution of the systems by using the method of Baker- Akhieser functions (see [3]). However, the 
resulting formulas appeared to be quite tedious, and it was not evident how to compare or identify 
them with the theta-function solution of Kotter. 
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Note that the latter was found in the classical manner, i.e., by a separation of variables and 
reduction of the equations of motion to quadratures, which have the form of the Abel-Jacobi map 
associated to a genus 2 hyperelliptic curve. The phase variables of the Kirchhoff equations have 
been expressed in terms of the separating variables in a quite symmetric but complicated way. Until 
recently, various attempts to check these expressions, as well as the reduction to quadratures made 
by Kotter, even using packages of modern computer algebra, were not successful. This even made 
some specialists to believe that the results of [16] are not reliable hence useless. 

One of the first step in verification of Kotters' calculations was made in [7J, where the Steklov- 
Lyapunov systems on e*(3), as well as their higher-dimensional generalizations, have been considered 
as Poisson reductions of certain Hamiltonian systems in a bigger phase space. The latter systems 
were shown to possess 2x2 matrix Lax representations in a generalized Gaudin form with a rational 
spectral parameter. This fact easily allowed to find separating variables, which coincided with those 
suggested by Kotter, and, as a byproduct, prove their commutativity with respect to the Lie-Poisson 
bracket on e*(3). A similar approach to the separation of variables was made in [22 . 

The main aim of the present paper is to reconstruct the rest of the results of the paper |16^1 . 

For our purposes we shall also use another set of phase variables which depend linearly on z,p. 
Namely, putting in (jSJ) successively s = b x , s = &2> s = 63 we obtain three independent quadratic 
integrals defining rank 3 quadrics in P 6 : 

(h - b 2 )(z 2 + b lP2 ) 2 + {bi - h){z 3 + b lP3 ) 2 = T{b x ) , 

(b 2 - h)( Zl + b 2P i) 2 + (62 - b 3 ){z 3 + b 2 p 3 ) 2 = T{b 2 ) , (10) 
(63 - b x ){z x + b 3Pl j 2 + (63 - b 2 )(z 2 + b 3P2 f = T{b 3 ) . 

Then it is natural to introduce new variables 

vi = \J (b 2 - b 3 )(bx - b 2 ) (z 2 + bip 2 ) 



v 2 = \J (b 2 - b 3 )(b 3 - bi) (z 3 + &ip 3 ) 



t'3 



= \J {b 3 - &i)(&i - b 2 ) (zi + b 2 px) 
va = \J (62 - b 3 ){b 3 - bi) (z 3 + b 2 p 3 ) 



v 5 = y/(b 3 - 61) (61 - 6 2 ) (zi + b 3 pi) 



v 6 = y/(b 2 - b 3 )(bi - b 2 ) (z 2 + b 3 p 2 ) 

which, in particular, imply 



(11) 



V 3 - V 5 Vl - V S V 2 - Vi 

Pi = -==—===, p 2 = 7= - = , p 3 = 



\/~S\/b 2 — 63 ' \f~S\/b 3 — bi ' y/S^/bi — b 2 ' 

S = (b 1 -b 2 )(b 2 -b 3 )(b 3 -bi). 

Then the integrals p0|) and (p, p) — Ji take the following compact form 

v 2 - v 2 = ^(h) / (b 2 - b 3 ) , 

vl-v^^/ibs-h), (12) 
v 2 - v 2 = V(6 3 ) / (61 - b 2 ) , 
(va - v 5 ) 2 (v x - v 6 ) 2 {v 2 - vi) 2 



b 2 -b z 63 - bi bi - b 2 



Ji(h - b 2 )(b 2 - b 3 ){b 3 - h) . 



The Steklov-Lyapunov systems written in terms of Vi, . . . , Vq , as well as the integrals (|12j) , are 
quite similar to those describing the reduction of the integrable geodesic flow on the group 50(4) 



1 Note that apart from the solutions of the Kirchhoff equations, Kotter also provided (although in an extremely 
brief form) the theta-solutions describing the motion of the group -E(3), that is, the components of the rotation matrix 
of the body and the trajectory of its center in space. We could not reconstruct these solutions. 
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with the diagonal metric II to the algebra so(4), which was considered in details in [UI2]- In fact, 
as was shown by several authors (see e.g., [5]), there is a linear isomorphism connecting the above 
system^. We shall use this property and the results of [5] to obtain theta function expressions for 
the sums and differences of Vi, which have an especially simple form. 



2 Separation of variables by F. Kotter. 

The explicit solution of the Steklov-Lyapunov systems in the generic case was given by Kotter in 
the brief communication |16j . where he presented the following scheme. 

Let us fix the constants of motion in ([9]), then the invariant polynomial (|8]) can be written as 

.F(s) = cn(s - ci)(s - c 2 )(s - c 3 ), c , ci, c 2 , c 3 = const. (13) 

Assume, without loss of generality, that b\ < b 2 < b 3 . Then one can show that for real z,p there are 
two possibilities: 

1) ci, C2, C3 are all real, then b\ < c\ < c 2 < c 3 < 63; 

2) ci is real and c 2 , C3 are complex conjugated, then b\ < c\ < 63 and either p — 5ftc 2 = 3?C3 < b\ 
or p > 63. 

Next, when no one of c a coincides with b±, 6 2 , 63, the level variety of the four first integrals of the 
problem (given by the coefficients at s 3 ,s 2 , s, s°) is a union of two-dimensional tori in R e — (z,p). We 
restrict ourselves to this generic situation, excluding the other cases, which correspond to periodic 
or asymptotic motions of the body. 

Let Ai, A 2 be the roots of the equation 

fiX)= J2 {ZjPk -_ Z b kPj)2 = 0, (ij,k) = (1,2,3), (14) 

i— 1 1 

where, when all c Q are real, 

AiG[6i,ci], A 2 e[c 3 ,6 3 ]. (15) 

Then for fixed Co, ci, c 2 , C3 the variables z,p can be expressed in terms of Ai, A 2 in such a way that 
for any s E C the following relation holds (see formula (7) in [TB]) 



3 



y/-(\i-c a )(\ 2 -c a ) ( yggOgAg y/gAgWAj 



( Cq - c ^)( Cq -c t ) I (Ai-6i)(A 2 -c a ) (A 2 -6»)(Ai-c a ) 
a— 1 ' 



where 



spi = Vcq — 5 = , (16) 

(Ai - A 2 ) ^ (Cq _ C(3)(Cq _ Ct) 



$(A) = (A-6i)(A-6 2 )(A-6 3 ), ^(A) = (A - Cl )(A - c 2 )(A - c 3 ) , (17) 



(i,j,fc) = (1,2,3), (a,/?, 7) -(1,2, 3). 

Setting in the above expression s — * 00 and s = 0, one obtains the corresponding formulas for Pi, Zi. 

Note that for real z^,^, in the case (1) (all c a are real), in view of the condition (fTS"]) all the 
expressions under the radicals in (116|) are non-negative. In the rest of the cases the roots can be 



2 On the other hand, one of the Steklov-Lyapunov systems on e*(3) can also be regarded as a limit of the system 
on so(4). 
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complex. For any a — 1,2,3, the branches of \f—{Ai — c Q )(A2 — c a ) in the numerator and the 
denominator of (fT6)) must be the same. 

Next, the evolution of Ai, A2 is described by the quadratures 

dAi , d\ 2 

di at , 



AidAi A 2 d\ 2 [ > 

— , H , = 02 at, 

R(X) = -$(A)^(A) 

with certain constants 81 , 82 depending on the choice of the Hamiltonian only. In other words, in 
the variables Ai,Aa the systems separate. 

Note that the paper [TB] does not describe explicitly how to find 61,62- They were calculated in 

mm- 

The above quadratures rewritten in the integral form 

r Xi dX r X2 d\ 

Ao 2/R(A) + Ao 2y/R(Xj~ Ul ' (2Q) 
XdX f Xl XdX 

U2 , 



Ao 2/R(A) 7a q 2/R(A) 
ui = Sit + u w , u 2 = S 2 t + U2o, (21) 

which represent the Abel-Jacobi map associated to the genus 2 hyperelliptic curve /i 2 = — $(A)^(A). 
Inverting the map (|20p and substituting symmetric functions of Ai, A2, /ii, /12 into (|16p . one finally 
finds z,p as functions of time. 

Everyone who had read paper |16j might be surprised by how Kotter managed to invent the 
intricate substitution (z,p) — * (Ai, A2, cq, ci, C2, C3) and to represent the result in the symmetric 
form (fl6|). Unfortunately, the author of [16] gave no explanations of his computations. Nevertheless, 
it is clear that behind the striking formulas there must be a certain geometric idea, which we try to 
reconstruct in the next section. 



3 A geometric background of Kotter's solution. 

Let (x\ : X2 : £3) be homogeneous coordinates in P 2 defined up to multiplication by the same 
non-zero factor. Consider a line I in P 2 = [x\ : X2 "■ X3) defined by equation 

ViXi + y 2 x 2 + 2/3^3 = 0. 

Following Pliicker (see e.g., [13]), the coefficients 2/1,1/2,2/3 can be regarded as homogeneous coordi- 
nates of a point in the dual projective space (P 2 ) . Now let h, 1% be two intersecting lines in P 2 with 

the Pliicker coordinates (y^, y^, 2/3 ), (2/1 ,2/2 :V^)- 

Then, for any constants A, n € C not vanishing simultaneously, the linear combination Xy^ + 
Hya^ are also Pliicker coordinates of a line l\ „ c P 2 . Hence, we arrive at an important geometric 
object, a pencil of lines in P 2 , i.e., a one-parameter family 1\ }I1 - It is remarkable that all the lines of 
a pencil intersect at the same point P £ P 2 called the focus of the pencil. 

Theorem 1. ([13 ) Let l\ tl± be a pencil of lines in P 2 defined by the Pliicker coordinates Xy^ +fJ-ya\ 
(A : fi) G P. Then the homogeneous coordinates of the focus are 

-r, / (1) (2) (1) (2) (1) (2) (1) (2) (1) (2) (1) (2)\ 

p = [V2 2/3 ~Va V2 • Vi Vs 2/3 Vi ■ Vi V\ ~ V2 Vi ) ■ 
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Next, consider the family of confocal quadrics in P 2 



Q(s) 



- b- 



= o 



(22) 



and a fixed point P = (X\ : X 2 : X3). Then one defines the spheroconical coordinates Ai, A2 of this 
point (with respect to Q(s)) as the roots of the equation 



0. 



A — bi A — b-2 A — 63 
Now, going back to the Steklov-Lyapunov systems, we make the following observation. 



Proposition 2. The separating variables Ai, A2 defined by formula \1J$ are spheroconical coordinates 
of the focus P of the pencil of lines in P 2 with the Pliicker coordinates z + sp = {z\-\- spi : z 2 + sp 2 : 
z 3 + s P3)j s G P with respect to the family of quadrics 



Proof. According to Theorem [TJ the homogeneous coordinates of the focus P are 

(z 2 p 3 - z 3 p 2 : z 3 px - zip 3 : zip 2 - z 2 pi) , (23) 

hence, the spheroconical coordinates of P with respect to the family (f2"2")l are precisely the roots of 
the equation (fT4)) , i.e., Ai,A2. □ 



Note also the following property: for a — 1, 2, 3, the line i a with the Pliicker coordinates z + c a p 
is tangent to the quadric Q a = Q(c a ). Indeed, setting in the right hand side of ([5]) s = c ai we obtain 



^(c a - bi)(zi + c a pif = , 



which represents the condition of tangency of the line l a and the quadric Q a . 

As a result, the following configuration holds: the three lines £i,£ 2 , £3 in P 2 intersect at the same 
point P and are tangent to the quadrics Qx, Q 2 , Q 3 respectively. An example of such a configuration 
is shown in Fig. 1. 

It follows that a solution z(t), p(t) defines a trajectory of the focus P on P 2 or on S 2 — {x\ +x 2 + 
x\ = 1}, and it natural to suppose that the Steklov-Lyapunov systems define dynamical systems on 
the sphere. Indeed, some of these systems were studied in [35] and were shown to be related to a 
generalization of the classical Neumann system with an additional quartic potential. 



Q 




Q, 



Figure 1: A configuration of tangent lines in R 2 
when the quadrics Q a are two ellipses and a hyperbola. 



||, || ) for the case 61 < a < b 2 < c 2 < c 3 < 63, 



G 



In the sequel our main goal will be to recover the variables z and p as functions of the sphe- 
roconical coordinates of the focus P, that is, to reconstruct the Kotter formula (fl6|) . Obviously, 
the solution is not unique: due to square roots in (fig)) , each pair (Xi, X 2 ), Xk 7^ b\,b2,b 3 gives 4 
points on P 2 , and for each point P that does not lie on any of the quadrics Q(c a ), 2 3 = 8 different 
configurations of tangent lines li,i 2 ,l 3 are possible (Fig. 1 shows just one of them). Thus, under 
the above generality conditions, a pair (Ai, A2) gives 32 different tangent configurations. 

Reconstruction of z,p in terms of the separating variables. Let (P 2 )* = (Gi : G2 ■ G 3 ) 

be the dual space to P 2 = (xi : X2 ■ x 3 ), (Gi being the Pliicker coordinates of lines in P 2 ). It is 
convenient to regard Gi also as Cartesian coordinates in the space (C 3 )* = (Gi, G2, G3). The pencil 
er(P) of lines in P 2 with the focus (|23[) is represented by a line in (P 2 )* or by plane 

7T = {(Z 2 P3 - Z 3 p 2 )Gl + (z 3 pi - Z 1 p 3 )G 2 + (Z\P2 - Z2Pl)G 3 = 0} C (C 3 )*. 

Consider the line tr(P) = {z + sp\s € K} C (C 3 )*. Obviously, {z + sp} C ir. Now let us use the 
condition for the three lines li,l 2 ,t-& defined by the points z + cip, z + C2P, z + c 3 p in (P 2 )* to be 
tangent to the quadrics Q(ci), Q(c 2 ), Q(c 3 ) respectively. Let V Q = (V a i,V a 2,V a3 ) C w, a — 1,2,3 
be some vectors in (C 3 )* representing these points, so that l a = {V a iXi + V a 2X2 + K*3^3 = 0}. Then 
we have 

z + cip - jtiiVi = , z + c 2 p - ^V 2 = , z + c 3 p- li 3 V 3 = (24) 

for some indefinite factors /i a - This system is equivalent to a homogeneous system of 9 scalar 
equations for 9 variables z a ,p a , Ha, a = 1,2,3. Thus the variables can be found up to multiplication 
by a common factor. Eliminating z,p from (|24|) . we obtain the following homogeneous system for 

Ml,M2,M3 

(c 2 - 03)14,1^1 + (c 3 - ci)V" a2 M2 + (ci - c 2 )V a3 fi 3 = , a = 1, 2, 3 , 

which has a nontrivial solution, since det ||V^,|| = (the vectors V a lie in the same hyperplane tt). 
It follows, for example, that 

Hi = flSi/ (c 2 - c 3 ) , fl 2 = l-iS 2 /(c 3 - ci) , H3 = ^^3/ (ci - c 2 ) , (25) 
£1 = V22V33 — V32V23 , £2 — V32V13 — V33V12 , ^3 = V12V23 — Vi 3 V22 , (26) 

fi =/= being an arbitrary factor. Substituting these expressions into (|24[) and using the obvious 
identity 

Si Vi + S2V2 + £ 3 V 3 = , 



after transformations we find 



(ci - c 2 )(c 2 - c 3 )(c 3 - ci) 



(ciSiVi + c 2 E 2 V 2 + c 3 i: 3 V3), (27) 



(c 2 c 3 S 1 V 1 + cic 3 S 2 V 2 + cic 2 X , 3V 3 ) . (28) 



(ci - c 2 )(c 2 - c 3 )(c 3 - ci) 



As a result, 



3 



z + sp = V" (c a s + cpc y )U a V a . (29) 

(ci -c 2 )(c 2 - c 3 )(c 3 - ci) ^ 

Now we express the components of V Q in terms of Ai, A 2 . Up to an arbitrary nonzero factor, 
they can be found from the system of equations 

3 

^1^+^2 + ^3X3=0, ^2(c a - bi)V*i = , a = 1,2, 3, (30) 

i=1 

which represent the conditions that the line £ a passes through the focus P = (x\ : x 2 : 2^3) and 
touches the quadric Q(c a ). 
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In the sequel we apply the normalization x\ + x\ + x\ — 1, which gives rise to expressions ([T%|) . 
For P ^ Q(c a ), this system possesses two different solutions, and for P £ Q(c a ) a single one (the 
line touches Q(c a ) at the point P). In the latter case we can just put 

V a i =Xi/ (c a - bi) . (31) 

Next, it is obvious that under reflection {x\ : x 2 : X3) — ► (— x\ : x 2 : X3), a solution (V a i : V a 2 : 
V a 3) transforms to (—V a i ■ V a 2 : V a3 ) (similarly, for the two other reflections). Let us seek solutions 
of equations (|30[) in the form of symmetric functions of the complex coordinates Ai , A2 such that 



1) for Xi = c a or A2 — c a (i.e., when P £ Q(c a )) there is a unique solution proportional to (|31[) : 

2) if Ai or A2 circles around the point A = c a on the complex plane A, the two solutions transform 
into each other; 

3) for Ai = bi or A2 = bi (i.e., when Xi = 0), V a i does not vanishes. 
Using the Jacobi identities 

{0, k<n-l 
I' k = n ~ 1 (32) 
Eoi, k = n, 

one can check that the following expressions satisfy equations (|30[) and the above three conditions 
_ / V<I>(Ai)(A 2 -c a ) VHX2)(Xi-c a ) \ _ V(A 1 -b i )(A 2 -6 < ) 



Then, using again the identities (|32p . we have 

(V a ,V^) = (A 2 -Ai)^ v /(A 2 -c Q )(A 2 -c^)- ^(A 1 - Cq )(A 1 - C(3 )^ 



(34) 



and, in particular, (V Q , V Q ) = (Ai — A2) 2 for a = 1, 2, 3. 

Next, substituting (I33p into ()26j) and applying the symbolic multiplication rule y/ab^/ac = a\/bc, 
we find the factors Z'q in form 

S a = (Ai - A 2 )xi (^-(Ai - c 7 )(A 2 - cp) - \J-(Xi - cp)(X 2 - c 7 )) , (35) 

(a, A 7) = (1,2,3). 
Further, putting (|33]l. (|35j) into (f29|). we obtain 

/i(Ai - A 2 )xi >r^, 

Zi + Spi = r- r- rXi ■ } {C a S + CttC^) 

ci-c 2 C2-C3 C3-C1) ;H 



^j»(Ai)V»(A 2 ) / / Ai-c T / Ai- C;3 \ A /$(A 2 )V>(Ai) / / A 2 -c T j Xi-c^ 
Ai— 6; \ V A 2 — e T y A 2 — I A 2 — 6,; I y Ai — c T y A 2 — c/3 



3 



x \ V^, ^ V-(A 1 -c Q )(A 2 -c et ) / yg(AQgAgj V$(A 2 )V-(Ai) \ , s 



which, up to multiplication by a common factor, coincides with the numerator in Kotter's formula 

To determine the factor \x in (|29p and in (|36p . we apply the condition (p,p) — cq which follows 
from (TT51) . Then, from (f2"T)) we get 



c Icx^xVi + 02^2 V 2 + c 3 r 3 V r |2 



,2 



(ci - c 2 ) 2 (c 2 - c 3 ) 2 (c 3 - Cl) 



(37) 
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Using the expressions (|34p . we obtain 



^C Q S Q V Q EE [4^a(Va,V a ) + 2c P CyEf,E 7 (V(l,Vj] 

a— 1 a— 1 

3 r / 

= (Ax - \ 2 fx\ 4(Ai - A 2 ) .^-(Ai - c 7 )(A 2 - C/J ) - \J—(\i - C/3 )(A 2 - c 7 ) 

+ 2c,gc 7 ^-(Ai - c 7 )(A 2 - ~Cfi) - yj ~{X\ - cp)(X 2 - c 7 )^j 
• ^-(Ai - c Q )(A 2 - c 7 ) ~ \J~(Xi - c 7 )(A 2 - c a )^j 

■ (V~( A 2 - C /3)( A 2 - c 7 ) - ^-(Ai - C /3 )(Ai - c 7 )^ . 

Simplifying the above expression and again using symbolic multiplication of square roots, one can 
verify that it is a full square: 



a=l 



= Xi(Ai - A 2 ) 4 [^{cp- c 7 )a/-(Ai - c Q )(A 2 - c a ) 



\a = l 



Hence, from ([57]) we find 



y/-{Xi - C Q )(A 2 - Ca) 



/' 



^ (Ca - Cp){c a - C 7 ) 



Combining this with (|36p . we finally arrive at (j!6[) . 

Thus, we derived the remarkable Kotter formula by making use of the geometric interpretation 
of the variables Ai, A 2 . We also note that the expressions (fT6|) are symmetric in Ai, A 2 . 



Remark 1. As noticed above, a disordered generic pair (Ai,A 2 ) gives 32 different configurations 
of tangent lines to the quadrics Q(ci), Q(c 2 ), Q{c%). Since the common factor fi in (|29|) is defined 
up to sign flip, we conclude that, according to the formula (116p . to each generic pair (Ai, A 2 ) there 
correspond 64 different points (z,p) on the invariant manifold (a union of 2- dimensional tori) defined 
by the constants Co, ci, 02,03. This ambiguity corresponds to different signs of the square roots in 
the Kotter formula. 

In the next section we shall use the expressions (fTBj) and the quadratures ([20]) to find explicit 
theta-functional solutions for the Steklov-Lyapunov systems. 



4 Explicit thet a- function solution of the Steklov-Lyapunov 

systems 

In order to give explicit theta-functions solution, we first recall some basic formulas describing 
inversion of the quadratures (fT^|) . We shall mainly follow the description given in 3, 4, 10J. Consider 
an even order hyperelliptic Riemann surface of genus g represented in the standard form 

r = {v 2 = (X-E 1 )---(x-E 2g+2 )} e c 2 (X,fi). 

In the sequel we shall regard T as its complex compactification obtained by gluing two infinite points 
oo_,oo + , where the coordinate A equals infinity. 

Consider also differential 1-form (differential) u = 0(r)dr on T, where r is a local parameter at 
a point Pef. A differential u> is called holomorphic if 4>(t) is a holomorphic function for any point 
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P. We choose the canonical basis of cycles Oi, . . . , a g , b\, . . . , b g on the surface T such that their 
intersections are of the form: 

en o dj ; = bi o bj ; = 0, Oi obj = Sij, i, j = 1, . . . , g, 

where 71 o 72 denotes the intersection index of the cycles 71 , 72 ■ 




Figure 2: 

An example of a canonical basis of cycles on T is shown on Figure 2. The parts of the cycles on 
the lower sheet are shown by dashed lines. 

Next, let u>i, . . . ,Q g be the conjugated basis of normalized holomorphic differentials on T such 
that 



1. 



The g x g matrix of 6-periods B^ = § b .0Ji is symmetric and has a negative definite real part. 
Consider the period lattice A = {2irjZ 9 + BZ 9 } of rank 2g in C 9 = (zi, . . . , z g ). The complex torus 
Jac(G) = C 9 /A° is called the Jacobi variety (Jacobian) of the curve G. 

Now consider a generic divisor of points P\ = (Ai,/ii), . . . , P g — (A 9 ,/U ff ) on it, and the Abel- 
Jacobi mapping with a basepoint Pq 



Pi 



Q + ■ 



(38) 



uj = lu>x, . 



,) T , z = ( Zl ,...,z g ) T eC 9 . 



Under the mapping, functions on S 9 T, i.e., symmetric functions of the points P\, . . . , P g are 2g-fold 
periodic functions of the complex variables Z\ , . . . , z g with the above period lattice A ( Abelian 
functions) . 

Explicit expressions of such functions can be obtained by means of theta-functions on the univer- 
sal covering C 9 = (zi, . . . , z g ) of the complex torus. Recall that customary Riemann's theta-function 
9{z\B) associated with the Riemann matrix B is defined by the seriea^l 



6{z\B)= exp({BM,M) + {M,z)), 
Mei, 9 

g 9 
(Af, z) = MtZi, (BM, M) = BijMiMj. 



(39) 



i=l 

Equation 9(z\B) = defines a codimension one subvariety Q £ Jac(r) (for g > 2 with singularities) 
called theta- divisor. 



3 The expression for 8(z) we use here is different from that chosen in a series of books on theta-functions by 
multiplication of 2 by a constant factor. 
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We shall also use theta-functions with characteristics a = [a\ , . . . , a g ) , (3 — . . . , g ), otj, 8j £ 
I, which are obtained from 8(z\B) by shifting the argument z and multiplying by an cxponcnQ 



a a 



01 ■ ■ ■ Pa 



(z) = exp{(Ba,a}/2+ (z + 2irj(3,a)} 6(z + 2irjf3 + Ba). 



Then for a pair of characteristics one has the following useful relations 



(z) = exp{(Ba', a') /2 + (z + 2ttj/3 + 2ttj/3', a')} I 



a + a 

All these functions possess the quadiperiodic property 



(z + 2ttj/3' + Ba'). (40) 



(*)■ 



(41) 



; (z + 2ttjK + BM) = exp(27rje) exp{-(BM, M)/2 - (M, z}}6 

e=(a,K)-(0,M), 

An important particular case is represented by theta-functions with half-integer characteristics 



V, 



and r ]tj =rn + r)j (mod I? 9 jl? 9 ) , A', A", rf{ e ^Z 9 /Z 9 



such that 



2ttj^ + B V ' 1 = / Co (mod A), 
2ttjA" + BA' = K (mod A), 



(42) 



K, € C 9 being the vector of the Riemann constants and Ei briefly denotes the branch point (^,0) 

on r. 



The half- integer characteristic ^ is odd (even) if 9 ^ (z) is odd (respectively, even). 

For the case g = 2 and for the chosen canonical basis of cycles 01,02, bi,b2 on T the above 
characteristics A, rji are 



A = 



V3 



1/2 


1/2 





1/2 





1/2 


1/2 






1/2 

0J ' 

1/2 

1/2 1/21' 



1/2 o x 
1/2 o y 



1/2 1/27 ' 



and, by convention, 776 is the zero theta-characteristic. Note also the property 

m + m + vs = m + m = A mod z 2g /z 29 . 



(43) 



(44) 



The six functions 9[A + rji](z), i = 1, . . . , 6 are odd, that is, 9[A + 77^] (0) = 0, whereas the other 
10 functions 9[A + rjij](z), i,jy^6 are even. In the case g = 2 no one of the latter functions vanishes 
at zero. 



The root functions. To obtain theta-functions solution for many problems linearized on Jaco- 
bians of hyperelliptic curves, one can apply some remarkable relations between roots of certain 
functions on symmetric products of such curves and quotients of theta-functions with half-integer 
characteristics, which are historically referred to as root functions. For the case of odd order hy- 
perelliptic curves such functions were obtained by Weierstrass and Rosenheim [331 HS]i see a ls° 

[3 CO- 

For our purposes it is sufficient to quote only several root functions for the particular case g = 2 
and the even-order hyperelliptic curve 

T = { M 2 = R{\)}, R(X) = (A — Ei) • • • (A — E e ). 

Let us introduce the polinomial U (A, s) = (s — Ai)(s — A2). 

4 Here and below we omit B in the theta-functional notation. 
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Proposition 3. Under the Abel-Jacobi mapping i38\) with g — 2 and the basepoint Pq = Eq the 
following relations hold 

U(X, Ei) s (A, - E,)(X 2 - E.) - - ^, ( ,f^g + , /a) . <«) 

Q = 2 j u, ^ = const, i — 1, ... ,6, 

JE 6 



\1-\2\ [Ei - \i){E 3 - Ai)(£ s - Ai) {Ei - \ 2 )(Ej - X 2 )(E S - A 2 ) y 
g[A + rh + 7?j + 7? a ](z) 6{A}(z - g/2) 8{A}(z + g/2) 
Kijs 6[A + T H ](z)9[A + T ]j ]{z)8[A + ri a \{z) 



(46) 



A1-A2 ^( J B i -A 1 )( J B i -A 1 )( J B s -A 1 ) {E i -X 2 )(E j - X 2 )(E S -A 2 ) / 

, 6>[A + r ? i + ? 7j +7 ?s ](z) 



(47) 



ys 0[A + r /s ](z) 

) «4' S = const > h s = 1, ■ ■ ■ , 6, i^ j ^ s^i, 

where, as above, t/q is the zero theta- characteristic and oo + ,oo_ are the infinite points of the com- 
pactified curve Y. The constant factors Kj, Kij s , K^- s depend on the moduli ofT only. 

Note that various expressions of symmetric functions of the A, /i-coordinates on an even hy- 
perelliptic curve were obtained in [11] on the basis of the Klein- Weierstrass realization of Abelian 
functions outlined in [3] and [5]. 

Sketch of proof of Proposition [3j The left and right hand sides of (|45p are meromorphic functions 
on Jac(r), which have the same zeros and poles with the same multiplicity. This implies that 
their quotient is an analytic function on a compact complex manifold without poles and therefore a 
constant. 

The root functions (|4^)) , (|4"?| can be deduced from the corresponding root functions for the case 
of odd-order hyperelliptic curve, by making a fractionally-linear transformation of A that sends the 
Weierstrass point E 2g+2 on Y to infinity. □ 

The constants Ki, K%j s , k'^ s can be calculated explicitly in terms of the coordinates E%, . . . , Eq and 
theta-constants by equating Ai, A2 to certain Ei and the argument z to the corresponding half-period 
in Jac(r) (see, e.g., [3]). 

Explicit solution. Now we are able to write explicit solution for the Steklov-Lyapunov systems 
by comparing the root functions (|45f . (J47J) with the Kotter expression (fl"6)) . 

Namely, let Y = |/i 2 = $ (A) (p (A)} where the polynomials <fi and ip are defined in (fTT|) and 
identify (without ordering) the sets 

{Ex, E 2 , E 3 , E 4 , E 5 , E 6 } = {b 1 ,b 2 ,b 3 ,a,c 2 ,c 3 }. 

By r]b { , i] Ca we denote the half-integer characteristics corresponding to the branch points (bi, 0), (c Q , 0) 
respectively, according to formula 



Theorem 4. For fixed constants of motion c\,c 2 ,c-$ the variables z,p can be expressed in terms of 
theta- functions of the curve Y in a such a way that for any s € C 

Zi + spi = —3 — l — — J - , {a, [3, 7) = (1,2,3), (48) 

where ki a , k 0a are certain constants depending on the moduli of Y only, and the components of the 
argument z are linear functions oft: 

z\ = CnSit + C\ 2 5 2 t + zio, z 2 = C 2 \b\t + C 22 S 2 t + zw, zw, z 2 q — const, C = A^ 1 (49) 
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A being is the matrix of a-periods of the differentials d\/ (m, XdX/ \i on T. 

Thus, we have recovered the theta- function solution of the systems obtained by Kotter in [TO] . 
The proof is given in the end of the section. 

Remark 2. In view of the definition of theta-function with characteristics, under the argument 
shift z — > z — JC the special characteristic A is killed and the solutions (|4"5)l are simplified to 

Zi + spi= —3 = — J , (a, 0, 7) = (1,2,3), (50) 

2^ a=1 k Oa 0[r] Ca \(z) 



where the constants ki a , ko a coincide with ki a , fcoa in (|48|) up to multiplication by a quartic root of 
unity. In each concrete case of position of bi, c a , one can also simplify the sums of characteristics in 
the numerator of (f50|) by using the relations (|44|) . 



Remark 3. In view of the quasi-periodic property (|41[) . when the complex argument z changes by a 
period vector in Jac(r), the theta-functions in (|48p . (|50[) are multiplied by generally different factors. 
Hence, the variables Zi,pi cannot be singlejvalued on the Jacobian variety T, and a simple accounting 
shows that they are meromorphic on Jac(r), the 16-fold unramified covering of it, obtained by 
doubling all the four period vectors in Jac(r). This implies that Jac(r) is also a principally polarized 
Abclian variety isomorphic to Jac(r). As follows from the structure of (|48|) . all Zi,pi have a common 
set of simple poles (the pole divisor), which we denote T> C Jac(r). 

The degree of the covering Jac(r) — > Jac(r) can also be found in another way: According 
to Remark 1, each generic pair (Ai,A2) corresponds to 64 different points (z,p) on the invariant 
manifold Jac(r). On the other hand, the same pair gives rise to 4 different points in Jac(r) defined 
by the divisors {(Ai, ±^/i? 6 (Ai) ), (A2, ±y/ Rq(X 2 ) ) }. Hence a generic point of Jac(r) corresponds 
to 64/4=16 points in Jac(r). 

Proof of Theorem^ The summands in the numerator of the Kotter solution (|16p . when divided by 
Ai — A2, can be written as 



V / -(Ai-c q )(A 2 - Cq ) / N /$(Al)^(A 2 ) V^WlAi 



(c a ~ cp)(c a - c 7 ) Ai - A 2 \(Ai - bi){X 2 - Ca) (A 2 - h)(Xi - c a ) 

_ s - c a \/-(Ai - cp)(\ 2 - cp)y/-(Xi - c 7 )(A 2 - c 7 ) 

(Cq - C/3)(c a - c 7 ) Ai - A 2 

Ml M2 



(Ai - 6»)(Ai - cp)(\i - c 7 ) (A 2 - 6 l )(A 2 - cp){X 2 - c 7 ) 



m = \/$(Ai)V'(Ai), /i 2 = \/ $ (^2)V'(A 2 ) . 

The right hand sides have the form of the root function (TfT)) . Hence, up to a constant factor, they 
are equal to 

J [A + T] C0 + T}^ + Tfo] (z) 



(« 

Next, in view of (1451), we obtain 



e[A + r, bi ]{z) 

e[A + r] bi ]{z) 



y/e[A}(z~q/2)9[A}(z + q/2) 



V-(A 1 -c a )(A 2 -c a ) = x a - rA1 / [A t^: r ] !!? (51) 

y/0[A](z-q/2)0[A](z + q/2) 



>Ci,>t a — const, 



\/-(Ai - c Q )(A 2 - c a ) _ ^ fc a^[A + ry Cc ,](z) 



^ (c Q - C/3 )(c Q -c 7 ) ^ A /0[A](z-g/2)fl[A](z + <//2)' 



(52) 
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Combining the above expressions, we rewrite the right hand side of (|1€>[) in the form 

9[A + i] bi ](z) ^ 3 k ia (s ~ c a ) 9 [A + r] Ct) + i] c ^ + i] bi ] (z) 



^9 [A] (z - g/2)9[A](z + g/2) ^ a ~ Y 0[A + 7fo](z) 

3 k Qa 9[A + Vca }(z) 



y/9[A}(z^q/2)e[A](z + q/2) 

which, after simplifications, gives 



Expressions (|49|) follow from the relation (toi,^) — C(dX/ /i, XdX/ /i), where, as above, uij are 
the normalized holomorphic differentials on T, which implies (u)i,W2) T = C(ui 1 U2) T , where U\,U2 
are the right hand sides of the quadratures (p?Tj) . □ 



5 The divisor of poles and the alternative form of the theta- 
function solution. 

The nice form of the Kotter solution (|4"5)) itself tells us a little about the structure of zeros and poles 
of Zi,Pi on the 2-dimensional Abelian variety Jac(T). It is possible however to give a quite detailed 
description of the set of common polesjaf these variables, called the divisor of poles T>. Obviously, 
V = {E«=i fc oa^[A + ffcJCO = 0} C Jac(r). 

Namely, for each a — 1,2,3, the zeros of 9[A + n Ca ](z) in Jac(r) form a translate <d a of the 
theta-divisor 9 by the half-period 2njri" a + 2Br/' Ca . Each translate passes via six half-periods, and 
61,82,03 have a unique common intersection in the origin (neutral point) O G Jac(r). This is 
depicted in Fig. 3 (a), where Q are shown as circles and the half-periods in Jac(r) as black dots. 
Hence, at z = O the denominator of (|4"5)) vanishes. Then, under the covering it : Jac(r) — > Jac(r), 
the preimage of O consists of all the 16 half-periods in Jac(r), which therefore belong to the divisor 
V. 

Note that translations in Jac(r) by a complete period V correspond to translation in Jac(r) by 
the half-period V/2. 

Now assume, as above, that bi < 62 < &3 and that (63,0) = Eg £ V is the basepoint of the Abel 
map (|38[) with g — 2. A further information about T> is given by 

Proposition 5. The divisor T> C Jac(r) is invariant under translations by the half-periods generated 
by 

Vi/2 = 2tt7< + 2B V ' bl , V2/2 = 2ttj< + 2B V ' b2 , K\ = n bl . (53) 

Proof. Choose a generic point q € T> and let z* be its projection onto Jac(r), which gives 

3 

f(z*) = Y / k Oa e[A + n Ca }(z*) = 0. 

a=l 

In view of the quasi-periodic property (|4ip and the half- integer characteristics f)43[) . under the 
translations z* — » z* + MVi + NV 2 , M, N e Z all the functions 9[A + rj Ca ] (z*) are multiplied by the 
same factor and therefore f(z* + MVi + NV 2 ) = 0. Hence the points z*/2 + MVi/2 + NV 2 /2 in 
Jac(r) also belong to V. 

One can also show that this does not hold for the translations by the other half-periods. □ 

Theorem 6. The denominator of the solution (j48[) admits the factorization 
3 

^2ko a e[A+r, Ca ](z) = cxp( X z+C)-0[A](z/2)0[A+ % ](z/2)0[A+? 7b2 ](z/2)0[A+7 ?bl +7 7b2 ](z/2) (54) 

ct=l 

with certain constants XiC- 
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The proof of the theorem is based on the fourth Riemann identity (see, e.g., [31 [TU]) and the 
theta- formulas of Frobenius and Thomae (see, e.g., [2lJ[l8]). Technically, it is quite tedious and for 
this reason we move it into Appendix. 

Now note that each of the 4 sets 

V Q = {6[A](z/2 \B) = 0}, V, = {9[A + %J(z/2 \B) = 0}, 
V 2 = {9[A + Vb2 ](z/2 \B) = 0}, V 3 = {9[A + Vbl + Vb2 ](z/2 \B) = 0} 

describe a translate of the theta-divisor, the genus 2 curve T embedded into Jac(r). Then, Theorem 
H] says that the pole divisor I? is a union of these translates, which are obtained from each other 
by shifts by the half-periods Vi/2, V2/2, and V3/2 = — Vi/2 — V2/2. The union passes through all 
the 16 half-periods in Jac(r). The action of the translations by Vi/2, V2/2, V3/2 in Jac(r) on the 
components (D ,T>i,T> 2 ,T> 3 ) gives respectively 

{V X ,V Q ,V 3 ,V 2 \ (2? 2 ,2?3,2?o,2?i), (V 3) V 2 ,V U V ). (55) 

All these properties are in complete correspondence with our previous observations about the 
divisor T>. 

Also, as was shown in [2] by applying the Kovalevskaya-Painleve analysis, the pole divisor with 
the same structure appears in the integrable flow on the algebra so(4) with the diagonal metric II, 
already mentioned in Introduction. This result of [2] about V equally holds for the Steklov-Lyapunov 
systems due to a linear isomorphism between them and the integrable flow on so(4). 

The intersection pattern for T> is shown in Fig. 3 (b), which we borrowed from [2J. Here the circles 
represent the translates T>j and the 16 black dots depict the half-periods. Under the projection it 
all the above half-periods are mapped onto O <G Jac(r). 




(a) (b) 



Figure 3: (a) Configuration of the translates & a in Jac(T). (b) The 4 translates of T in Jac(T) forming the pole 
divisor T>. 



Solutions for the variables Vk- Let us choose the origin of Jac(r) at one of the four triple intersec- 
tions of T>j and denote for brevity the four theta-functions in ([54]) as 9q(z/2), 6\{z/2), 02(z/2), 9 3 {z/2). 

Now we show that theta- function solutions for the new phase variables v\ , . . . , vq introduced in 
(jTTJ) have a rather specific and compact form. Namely, as follows from expressions (j48|) and (fTl]) , the 
functions v\ + v 2 and v\ — v 2 may have only simple poles at most along the components of the divisor 
V. On the other hand, the form of the integrals (fi"2|) imply the following remarkable property: the 
poles (the zeros) of Vi + v 2 are the zeros (resp. the poles) of V\ — v%. Since both functions are 
meromorphic on Jac(r), none of then can have simple poles along only one component Dj. This 
necessarily implies that v\ + v 2 has poles along two certain components T>j 1 , T>j 2 and zeros along the 
other two components T>j 3 , 2? J4 , and vice versa for v\ — v 2 ■ 
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The same observations hold for the pairs (v 3 + U4, v 3 — V4) and (^5 + vq, v$ — vq). Note also that 
functions from different pairs cannot have the same poles, since in that case they would also have 
the same zeros and their quotient would be constant, which is not true. 

Now let us fix the origin of Jac(T) at one specific triple intersection of T>j such that the 3 functions 
v\ + v 2 ,v 3 + V4,, U5 + vq have a common pole along the component T>q. In this case the following 
proposition holds. 

Proposition 7. The theta- function solutions for the phase variables Vk have the form 

0i(z/2)fl 2 (z/2) _ _ fl (z/2) 6 3 {z/2) 

Vl 2 Xl e (z/2)9 3 (z/2)' Vl V2 X2 6 x {z/2)e 2 {z/2y 
9 2 (z/2)6 3 (z/2) g (z/2)g!(z/2) 
9o{z/2)0-l{z/2) 02{z/2)e 3 (z/2) 
_ gi(g/2)^(g/2) _ _ 8 (z/2)8 2 (z/2) 

V5+V6 - X5 e (z/2)6 2 (z/2Y V5 V6_X6 6 x {z/2)6 3 {z/2y 

^(61) tM Mh) (VA 

Xi,X3,X5 = const, x 2 = 77 j—, — , X4 = 77 — , *6 = 77 r"x — > ( 57 ) 

(O2-O3JXI (O3-O1JX3 (Ol-0 2 )X5 

where z — [z\, z 2 ) depend on t according to (|49D . 

Proof. First, note that the functions (|5rj|) have the same structure of zeros and poles, as prescribed 
above. Next, as follows from the Kotter formula (|16p and theta-solutions (|4"5)) . the translations by 
the period vectors Vi, V 2 , Vi + V 2 in Jac(r) generate the involutions 



cri 

03 = o 2 o 0\ 



(z 1 ,p 1 ,Z 2 ,p 2 ,Z 3 ,p 3 ) H-> (zi,Pl,-^ 2 ,-P2, 23,^3), 

(zi,pi,z 2 ,p 2 ,z 3 ,p 3 ) 1 ► (-zi,-pi,-z 2 ,-_p 2 ,z 3 ,p 3 ), 
(zi,pi,z 2 ,p 2 ,z 3 ,p 3 ) 1 ► (-zi,-pi,z 2 ,p 2 ,z 3 ,p 3 ), 



which, in view of (jTTJ) , gives rise to the transformations 



CI 


v 2 


+ Vl 


0-2 


v 2 


+ «1 


0- 3 


v 2 


±«1 



V 2 -Vl, V 4 ±V 3 < >U 4 ±W 3 , W5+«6< >V 5 -V 6 , 

V 2 -Vl, V4+V 3 < >V A -V 3 , V 5 ±V 6 < >-(v 5 ±V 6 ) 

V 2 ±Vx, V 4 +V 3 < >V4-V 3 , v 6 +v 5 < >v 6 -v 5 . 



Now observe that the relations ([55]) are invariant under the action of Oi on the left-hand sides and 
the corresponding transformation of 9o(z/2), . . . , 8 3 (z/2) under the action ([55]) . Moreover, one can 
check that the left- and right hand sides of (|56[) are multiplied by the same factors under the shift 
of z by any period vector of Jac(r). This proves |56|) . 

The relations l[57]) between the constants \i follow from the first 3 integrals in (fT2"|) . □ 

The constants xi> X2> X3 can De calculated explicitly in terms of bi, c a and theta-constants of T. 
As follows from the solutions (|56|) . the product («i + u 2 )(i>3 + W4) and the other two similar 
products have double poles along T)q only: 

/ w v ^2(^/2) 

(Wl+V2)(V3+V4)=52 ^72)' 

»l(*/2) , v/ s Sl{z/2) 

(V 3 + «0(t* + «,) = 53^, («! + V 2 )(V 5 + V 6 ) = 9l ^ y 

9i,92,g 3 = const. 

Analogs of some of these expressions were obtained in paper [5] in relation with separation of variables 
for the integrable system on so(4) with the diagonal metric II. Due to the linear isomorphism between 
this system and the Steklov-Lyapunov systems, the separating variables presented in [9] can also be 
regarded as new separating variables for (j4]), (J5j) . 
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6 Conclusive Remarks 



In given paper we gave a justification of the separation of variables and the theta-function solution 
of the Steklov-Lyapunov systems obtained by F. Kotter [IB]. Using the results of dJ[5], we also 
presented such a solution for an alternative set of variables, which have a simpler form. 

On the other hand, there exist several nontrivial integrable generalizations of the systems: the 
first of them was discovered by V. Rubanovsky [TI5] and describes a motion of a gyrostat in an ideal 
fluid under the action of the Archimedes torque, which arises when the barycenter of the gyrostat 
does not coincide with its volume center. In this generalization the Hamiltonian of the Kirchhoff 
equations, apart form quadratic terms, contains linear (gyroscopic) terms in K,p. Under the change 
of variables ([3]), the gyroscopic generalizations of the systems (HJ), ([3]) take the form 

i = z x (Bz — g) — Bp x (Bz — g) , p — p x (Bz — g) 

and, respectively, 

z = p x (Bz — g) , p = p x (z — Bp) , 

where g — (g%, g 2 , <?3 ) T is an arbitrary constant vector related to the angular momentum of the rotor 
inside the gyrostat. 

Following [12] . these systems admit the following generalizations of Kotter's Lax pair with an 
elliptic spectral parameter 

L(s) = [L(s),A(s)}, L(s),A(s) e so(3), seC, 
L(s) a p = e a p^U/s - by (z 7 + sp 7 ) + g 7 / yfs - b^j , 

^-(s) Q /3 = £ <xPi-\j (s - b a ){s - bp) (b y z^ - p 7 ) , resp. A(s) a/3 = e aj3l \J [s - b a )(s - bp) p 7 , 

which provides a sufficient set of constants of motion and makes possible to obtain theta-function 
solutions. Like in the case of the Steklov-Lyapunov systems, generic invariant manifolds of the 
Rubanovsky systems are two-dimensional tori, which can be extended to afhne parts of Abelian 
varieties. However, as we plan to show in a forthcoming publication, an explicit integration of the 
latter systems appears to be more complicated, and the Abelian varieties are not Jacobians of genus 
2 hyperelliptic curves, but Prym subvarieties. 

The problem of separation of variables for the Rubanovsky systems is still unsolved. 
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Appendix. Proof of Theorem M- 

The proof is based on the fourth Riemann identity (see, e.g., [3l 110)) 



%l)% 2 )%3)%4) = 2E <? 



(toi)i 



{w 2 )0 



{w 3 )t 



(Wi), 



(58) 



where the summation is over all the half-period characteristics 
(in our case g = 2) are related as follows 



and the arguments yj , wj S C 9 



(wi w 2 w 3 w 4 ) = (yi 2/2 2/3 2/4) T, T = i 



/ll 1 1\ 

1 1 -1-1 
1-11-1 

\i -1-1 1/ 
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Up to multiplication by a simple exponent of z, the theta- product in (|54[) can be written as 

6{z'/2) 8(z'/2 + Vi) 6(z'/2 + V 2 ) 6(z'/2 + Vi + V a ), (59) 
where z' = z + 2/C, i.e., the translation by the complete period in Jac(T), and Vi,2 are the periods 



defined by (j53|) . In view of the identity (|58|) . the product ([59]) gives the following sum of 16 theta- 
products: 



E 

2(a,/3)e(Z 2 )4 



z' + 



Vi + v 2 





a 




a 




a 


)' 


A 




A 




A 



(0). 



(Note that in each product the variable z' enters only once.) 

Next, in view of the property (|40[) . this sum can be written as product of an exponent of z and 
the sum 



E 

2(a,/3)S(Z 2 ) 4 L 



a + a 
P + (3' 



CO 



a + a 

[3 + [3> 



Vi 
2 



a + a 
(3 + 13' 



2 



a + a 
(3 + (3' 



V 1 + V 2 



2ttj f3' + Ba' = (Vi + V 2 )/2, 
which, under the corresponding re- indexation, reads 



E £q ^ 6 ' 

2(a,/3)G(Z 2 )4 
6 



Vi + v 2 



~ 0[A + Vij ](z') 8[A + Tm] (?f) 8[A + mj ] C^f) 6[A + m 



l<i<j<5 



2 

Vi+V 2 
2 

V! + V 2 



(60) 



where e. 



a,/3 



-1 if 8[a 0\(z) is odd and +1 otherwise, and, as above, rjij = rjj + r]j mod Z 2 /2 



In fact, most of the theta-constants in (f6T)|) are proportional to 8[A + T]i](0), i — 1, ... ,6 and 
therefore vanish. Namely, in the first sum in the right hand side of (|60[) all the theta-constants 
are non-zero if and only if rji is different from T] bl ,ri b2 , and 776 = 0. In the second sum, if rji or r]j 
coincides with rj bl or rj b . 21 then either the first or the second theta-constant is zero. Otherwise, if 
{ViiVj} n {^1,^2} = $i then, in view of the relations (|44|) . the third theta-constant is proportional 
to 8[A + ?/fc](0), for a certain k e {1, . . . , 6} and, therefore, equals zero. 

Since for the case of genus 2 there are no even theta-functions which vanish for zero value of the 
argument (see [3]), one concludes that the above sum contains only 3 non-zero theta-products: 



Y, 6[A + Vi ]{z')e[A + m ] 

3 



8[A + rji] 



Y 0[A + ffcJCO 6[A + Vca ] -L \ 9 [A + r,, 



8[A + m ] 



V1+V2 
2 

Vi+V 2 



Now, assume (for the moment) the following ordering of the Weierstrass points: 
£1=61, E 2 =ci, E 3 = b 2 , Ei = c 2 , E 5 = c 3 . 



(61) 



Then, in view of (|40|) and the identities (|44[) . the above sum, up to a constant common factor, can 
be written as 



e^[A + ry cl ](z')^[ 7?6l + 7?C2 ](0)^[7 ?fe +77 C2 ](0)e[j 7ci +7 ?C3 ](Q) 
+e 2 8[A + r) C2 ](z') 8[ Vbl + r) Cl }(0) 8[ m2 + Vci }(0) 8{r, C2 + Vc3 }(0) 
+e 3 9[A + r, C3 }(z')e[r) b2 }(0)6[ Vbl ](Q)e[Q}(0), 



(62) 
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where now e a are certain quartic roots of 1. 

Now we are going to show that the denominator in the theta-function solution (|48|) coincides 
with (|62|) up to multiplication by an exponent of z'. 

Namely, in view of (|52[) , the sum Y2a=i ^0a9[A + r) Ca ](z) can be written as a product of 

const y/0[A](z - q/2)9[A](z + q/2) 

and the expression 

g _ \/-(Ai - ci)(A 2 - ci) / c 2 - c 3 | v/-(Ai - c 2 )(A 2 - c 2 ) /c 3 - c x { ^/-(Ai - c 3 )(A 2 - c 3 ) 



a/( c i _ c 2)( c i = c 2) V ci - c 2 ^/(ca - ci)(c 2 - c 3 ) V ci - c 2 ^(cs ~ ci)(c 3 - c 2 ) 

(there is no second radical in the third summand !). Now, make the projective transformation 
A — > 1/ = A/(A — _Eg) = A/(A — 63), which sends the Weierstrass points c a , b%, 6 2 , 6 3 on F to c a , b\, 6 2 , 
and 00. The two infinnite points over A = 00 are mapped to 2 points over 1/ = 1. This change leaves 
the sum Q almost invariant: it becomes the product of const/ \J (y\ — l)(f 2 — 1) and the sum 



q _ - ci){v 2 - ci) / c 2 - c 3 | c 2 )(i/ 2 - c 2 ) /c 3 - ci ( - c 3 )(z/ 2 - c 3 ) 



- c 2 )(ci - c 2 ) V ci - c 2 ^/(ca - ci)(c 2 - c 3 ) V ci - c 2 y/ (c 3 - ci)(c 3 - c 2 ) 



Under the Abel map (J38J) , the radicals in Q can be expressed completely in terms of the theta- 
functions and theta-constants of F: Applying the theta-formulae of Frobenius and Thomac for the 
case when one of the Weierstrass points of the curve lies at infinity (see, e.g., [2UHE]) anc ^ keeping 
the ordering (|6ip . we have 



y/ H-Cq )(v2~c a ) g A + 9 A + Vc J (z 

y(c~r^)^r-r7ry 0[A + r? cl + ?7 ca +7? C3 ](0) 0[A](*) 

c 2 - c 3 = ^ 9[q C2 +7? fe2 ](0)6>[77 C2 +7? fcl ](0) ^ 

ci - c 2 % Cl + ?7 C2 + r; C3 + ?7h 2 ](0) 9[rj Cl + ?y C2 + ry C3 + %J(0) ' 
C3-C1 % Cl + 7/ 5l ](0) % Cl + 77 h2 ](0) 

— 03" 



ci - c 2 ' 6[q Cl + r] C2 + r] C3 + r?6 2 ](0) 9[i] Cl + T) C2 + rj C3 + %J(0) ' 

where r\ Ca , rj bl , rj b2 are the same as above and Qi are the appropriate quartic roots of 1 . Lastly, we 
have 

r — ygjAj (z-q/2)9[A}(z + q/2) 

a/(i/i - l)(i*2 - 1) = const 9[A]{z) ' ^ ' 

where q is the same as in (|45[) . or, in terms of the new coordinate v on T, q/2 = J 



(1,0) 



Combining the above expressions, we see that in the quotient Q = Q / \J (y\ — \)[y- 2 — 1) the 
term 0[A](z) is canceled and in the product Q\j9\A\(z — q/2)9[A](z + q/2) the square root (|64|) is 
canceled. Now, simplifying the theta-characteristics in (|63p by using (|44[) and ignoring common 
constant factors, we eventually find 



3 

const fc a^[A + ?7 c J(z) 

a=l 

= £!0[A + Vci ](z) 0[ Vci + Vc3 }(0) 0[ Vc2 + Vbl }(0) e[ Vc2 + n b2 ](0) 
+ E 2 9[A + i lc2 ]{z)9[r jc2 +J]c 3 }(0)e[r ]ci + ^(0) 0[t? Ci + ^ 2 ](0) 
+ e 3 9[A + Vc3 ](z) 9[ Vb2 }(0) 9[r, bl }(0) 9(0), (65) 

Ej also being certain quartic roots of 1. The latter expression have the same structure as the sum 
((62)) . Lastly, note that under the shift of z by an appropriate complete period in Jac(F) the roots 
£i can be made proportional to any combination of roots e a in (|62p . (This corresponds to choosing 
an appropriate origin in Jac(F).) Hence, we proved the theorem for the chosen ordering (|6ip . 

To complete the proof for the other possible orderings of bi , c a it remains to modify the theta- 
characteristics in (|62|) . (|65"|) . □ 
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